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Abstract 



We present simulation data for the motion of a polymer chain through 
a regular lattice of impenetrable obstacles (Evans-Edwards model). Chain 
lengths range from N = 20 to N = 640, and time up to 10 7 Monte Carlo 
steps. For N > 160 we for the central segment find clear i 1/,4 -behavior as an 
intermediate asymptote. The also expected i 1,/2 -range is not yet developed. 
For the end segment also the i 1//4 -behavior is not reached. All these data 
compare well to our recent analytical evaluation of the reptation model, which 
shows that for shorter times (t < 10 4 ) the discreteness of the elementary 
motion cannot be neglected, whereas for longer times and short chains (N < 
100) tube renewal plays an essential role also for the central segment. Due 
to the very broad crossover behavior both the diffusion coefficient and the 
reptation time within the range of our simulation do not reach the asymptotic 
power laws predicted by reptation theory. We present results for the center- 
of-mass motion, showing the expected intermediate t 1//2 -behavior, but again 
only for very long chains. In addition we show results for the motion of the 
central segment relative to the center of mass, where in some intermediate 
range we see the expected increase of the effective power beyond the i 1//4 -law, 
before saturation sets in. Analysis and simulations agree on defining a new 
set of criteria as characteristic for reptation of finite chains. 
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I. INTRODUCTION 

An understanding of the motion of a chain molecule in a surrounding of impenetrable 
obstacles is of great interest in the physics of polymer melts or dense solutions as well as 
for polymers diffusing through gels. With special regard to the latter system De Gennes 
suggested the reptation model |1J] . Basic to this model is the observation that the crosslinked 
structure of the gel for short times restricts the motion of the macromolecule to a tube defined 
by its initial configuration. The motion proceeds by curvilinear diffusion of little wiggles of 
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'spared length' along the tube. The destruction of the initial tube ('tube renewal') is due 
to the motion of the chain ends. These may draw back into the tube, thus shortening the 
tube and creating new wiggles of spared length, or they may unfold and thus destroy spared 
length, thus prolonging the tube in some random direction. This is the natural thermal 
motion of a flexible chain between topological constraints. 

Proposed originally for motion through rigid gels, this model extensively has been applied 
to melts or dense solutions M. It is generally accepted as a basic scenario of polymer 
dynamics. A critical examination ||, however, shows that the experimental or computer- 
experimental evidence for the quantitative reliability of the model is not particularly strong. 
Searching for the asymptotic power laws predicted by reptation theory one typically finds 
a range of exponents differing from the predictions, a finding often interpreted as crossover 
behavior from Rouse-type motion to reptation. This is little more than an excuse, since 
-, with the exception of Doi's theory of the melt viscosity j| -, no effort seems to have 
been spent to really work out the predictions of the reptation model beyond (intermediate) 
asymptotics. Thus, since not only the experiments but also the simulations mostly are 
concerned with melts or with an immobile, but disordered configuration of obstacles, it is 
not at all clear whether the results reflect intrinsic properties of reptation or are dominated by 
other mechanisms like entropic traps in disordered systems or relaxation of the surrounding 
in melts. 

To proceed we need precise knowledge of the quantitative implications of the repta- 
tion model in the (computer-) experimental range of time and chain length. We therefore 
analytically have worked out detailed quantitative predictions of the model, and we have 
carried through extensive simulations. All our work is concerned with the original reptation 
scenario: motion of a discrete chain through an ordered lattice of impenetrable obstacles, 
which confines the internal motion of the chain to a very narrow tube. Our analytical work 
concentrates on the motion of individual beads. As long as the bead does not leave the 
original tube, its motion can be calculated rigorously. Tube renewal and thus the motion 
of the chain ends can be treated only approximately, and we use an approximation inspired 
by random walk theory. Details and results of our analytical work may be found in the 
preceeding paper 0. 

The present paper is devoted to our simulations. In Sect. II we introduce a Monte Carlo 
model, which first has been proposed by Evans and Edwards ||. We also briefly describe 
our analytical model as far as needed for some of the arguments to follow, and we discuss 
the relation among the Monte-Carlo and analytical models. In Sect. Ill we review results of 
previous simulations of the Evans-Edwards model, compare to corresponding results of our 
simulations, and discuss the relevant time scales. Sects. IV and V are devoted to a detailed 
comparison with our theory, where in Sect. IV we treat the motion of the central bead 
inside the tube, and in Sect. V we are concerned with tube renewal. Quantities involving 
the center-of-mass, for which at present we have no new analytical results, are discussed in 
Sect. VI. Sect. VII summarizes our findings. Preliminary results of our analytical work and 
our simulations have been published in 

II. MODELS 
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A. Monte Carlo model 



The Evans- Edwards model || considers the chain configuration as a random walk of 
jy\MC) _ i s t e p S ('segments') on a cubic lattice. The lattice constant £q henceforth is taken as 
the unit of length: £ — 1. The configuration is fixed by giving the positions {r 1; . . . , r N ( M c) } 
of the endpoints of all segments ('beads'). The length of segment j equals £ , by construction, 
|r 3 - + i — Tj\ = £q. The obstacles form a second cubic lattice, of lattice constant m ■ £ , 
placed such that its lattice points coincide with centers of the cells of the first lattice. 
The edges of this lattice are considered as impenetrable. Ref. uses m = 1, 2, . . . , 10, 
but we consider only m = 1, thus taking the tube as narrow as possible. This should 
show reptation in clearest form. As illustrated in Fig. 1, this model eliminates all kink-type 
motions of the chain and leaves only hairpins, i.e., subsequent segments of opposite direction: 
rj + i — Yj = — (r? — free to move. Of course also the end beads can move freely. 

Clearly with regard to the static properties this model is identical to a simple random- 
walk chain. With our choice of the narrowest tube, m — 1, also the dynamics is most 
simple. The obstacle lattice comes into play only implicitly in restricting the motion to that 
of hairpins and chain ends. We start with an initial random walk configuration of the chain. 
In one elementary step we randomly choose one bead. If it happens to be the tip of a hairpin 
or a chain end, it is moved with probability 1/6 to one of its 6 possible positions (including its 
original position, of course). This completes the elementary move. Monte Carlo time t (MC ) 
is measured in (on the average) one attempted move per bead. The simulations extended 
to t (MC 1 = 10 8 , and chains of lengths N^ MC ^ = 20,40,80,160,320,640 were used. We 
measured in each run correlations over time intervals t[ MC ^ — tQ MC ^ < 10 7 , averaging over 
^(mc) (> mov j n g average'). In addition the data were averaged over up to 40 independent 
runs. This is important in particular for the longer chains, where the equilibration time 
T^ MC ^ of the hairpins comes close to the total time of the run (see Sect. 3.2). 

For the longest chains (N = 320, 640) and largest times (t( MC ^ « 10 7 ) the standard 
deviation of our data reaches 6 %. Due to the moving time average it rapidly decreases with 
decreasing chain length and time, being less than 3 % for t^ MC ^ < 10 5 , for all chain lengths. 



B. Analytical model 

We use a version of De Gennes' reptation model [0,0, discrete in time and space. The 
tube is taken as a chain of N segments, connecting beads numbered 0, 1, . . .N. Particles, 
each representing a spared length £ s , are sitting on the beads of that chain. The average 
density of these particles is po- The particles randomly and independently hop among 
neighbouring beads, with hopping probability p. They do not interact, so that a given 
particle does not feel the presence of the others. If a particle moves over a bead, it drags it 
along and displaces it by a distance £ s along the tube. (For hairpin motion this is illustrated 
in Fig. 1 b of 0.) Reaching a chain end (bead or N), a particle is absorbed by a virtual 
reservoir. These reservoirs also randomly emit particles at a rate adjusted to pertain the 
average density po- They serve to ease the analysis of the in principle grand canonical 
problem. 

To establish the connection to the physical motion of the beads we note that the dis- 
placement along the tube of some bead j within time interval t is given by £ s \n(j, t)\, where 
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n(j,t) gives the number of particles having passed bead j from one direction, subtracted by 
the number of particles coming from the other direction. Since the tube conformation itself 
is a random walk in space, bead j in space has moved an average distance 

g 1 (j,N t t) = ((r j (t)-T j (0))*) 

= £sHJM ■ (2-1) 

(Cf. Eq. (I 2.8); in the sequel ref. |J will be refered to as I.) In Eq. (2.1) the pointed 
brackets denote the average over the chain configurations, and the bar indicates the average 
over particle diffusion. 

Eq. (2.1) holds as long as bead j stays in the initial tube. Tube renewal is driven by 
the emission and absorption of particles by the reservoir. Emission of a particle shortens 
the tube by i s at the end considered. Thus within time interval t the tube from end zero is 
destroyed up to the bead j < = £ s n max (t), where (— n max (t)) is the largest negative fluctuation 
in the occupation number of reservoir within time interval t. Since at time t the tube on 
the average has been shortened by £ s n max (t) steps and then rebuilt by another £ s n max (t) 
randomly chosen steps, we find for the motion of the endsegment (cf. Eq. (I 2.12)) 

Pi(0, N,t) = 2£ s (t) , (2.2) 

valid for times t smaller than the tube renewal time T3. Combining these considerations 
we find the somewhat complicated expression (I 2.13) for the motion of an arbitrary bead 
including tube renewal effects. 

All these are exact expressions within the frame of our model, holding as long as the orig- 
inal tube is not destroyed completely. It turns out that Eq. (2.1) can be evaluated rigorously, 
whereas Eq. (2.2) as well as the tube renewal effects on the motion of an arbitrary bead 
can be handled only approximately. The stochastic process n(0,t), giving the occupation of 
reservoir 0, is correlated by the fact that a particle emitted may be reabsorbed at some later 
time, the decay time of the correlation being given by the time T 2 a particle needs to diffuse 
over the whole chain. This correlation renders an exact evaluation of n max (t) impossible. As 
explained in detail in I, sect. V we evaluate n max (t) in a 'mean hopping rate' approximation, 
calculating the contribution to n max {t) of a time step s, < s < t as the contribution of 
an uncorrelated process with properly adjusted hopping rate. In a similar spirit we have 
constructed an approximation for the tube renewal effects on arbitrary beads. Our explicit 
expressions for gi(j,N,t) will be recalled later in the context of data analysis. 

C. Relation among the models 

Obviously the particles of the theoretical model roughly correspond to the hairpins, 
and for long chains and large time, where hopefully the influence of the microstructure is 
negligible, we expect both models to yield identical results. In practice, however, it is not 
clear whether the experiments reach such a universal regime, and a more detailed discussion 
of the relation among the models is appropriate. 

We first consider the chain lengths. The endbeads of the MC-chain correspond to the 
particle reservoirs, and a hairpin absorbs two beads. A hairpin thus effectively walks along 



4 



a chain of N^ MC > — 4 interior beads. This must be compared to the theoretical model, where 
a particle hops along a chain of iV — 1 interior beads. Thus we should identify 

N = N {MC) - 3 . (2.3) 

For the shorter chains (N( MC ^ < 100), this correction cannot be neglected. 

Identifying hairpins and particles we should take the spared length £ s = 2. The density 
p is less well defined. For the simple random-walk type MC-chain it is not hard to determine 
the full statistics of the side-branches, i.e., of tree-like structures in which each lattice bond 
is occupied by an even number of segments. In the limit of long chains the average total 
number of segments in such side branches amounts to [0] N^ MC ' /3, whereas the average 
number of simple hairpins (of two segments each) tends to N^ MC ' /9. Thus the remaining 
segments are contained in larger side-branches, which can be seen as a result of 
a fusion of simple hairpins. In the particle picture this corresponds to an interaction of 
particles sitting on the same bead, and with this interpretation we should choose p a ~ 1/6. 
However, also other complications must be noted: For a hairpin lying on the chain like in 
the right part of Fig. 1, the separation of the configuration into hairpin and backbone of 
the chain is not unique, this configuration in fact showing two mobile points. Thus these 
considerations suggest an order of magnitude for £ s ,p Q , rather than giving precise values. 
Taking i s , p Q as fit parameters, we in Sect. 4 will find that the data rather precisely determine 
the combination 

f s Po = 1.23 (2.4) 

which, as shown in I, sect. IV is the only combination relevant in the universal large time 
regime. With E 2 s p fixed, a range of values 0.15 ^ p a ^ 0.3 yields equivalent fits, and by 
convention explained below we choose 

Po = 0.22 , (2.5) 

leading to 

£ s = 2.364 . (2.6) 

These parameters are of the expected order of magnitude, but they also show that the iden- 
tification of particles and 'free' hairpins should not be taken to literally. Rather the hairpin 
motion is renormalized by interaction effects, the particles representing 'quasi-hairpins'. 

We finally consider the relation among the time scales. The theoretical results, considered 
as function of 

t = pt , (2.7) 

for p ■ t ^ 1 essentially are independent of the hopping rate p. We by convention take 
p = 1/5, and we henceforth always will use the variable i. The relation among i and t^ MC ' 
defines the time scale r: 

t = rt^ . (2.8) 
A fit to experiment (see Sect. 4) fairly precisely fixes r at a value 
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r = 6.092 • 10" 2 



(2.9) 



Thus about 17 MC moves correspond to the displacement of a particle by one step. This 
again is a reasonable result, since following the microscopic motion of a hairpin we may 
estimate that on the average of the order of 10 moves are needed for a hairpin to jump from 
one segment to the next. 

Having discussed the relation among the parameters of the theoretical and the MC-model 
we still need to consider the measured quantity 

{MO^^MO^MC^ = /(j^iMC)) _ r . (Q)) 
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Let j be some interior bead of the MC-chain. With probability pn — 1/9 it sits in the tip 
of a simple hairpin, a configuration which in the theoretical model effectively is projected 
down to the base of the hairpin. Taking into account only simple hairpins we thus find for 
the relation of g[ MC ^ to the g\ of the analytical model 

g[ MC \j,N^ c \t^) = {l- PH f 9l {j - 1,NW-S,t) 

+ 2p H (l - p H ) \ 9l (j - 1, - 3, t) + 1 



Ph 



where we took into account the relation among N^ MC "> and N as well as the different counting 
of the beads. This shows that g[ MC ^ and g\ differ by an additive contribution cq 

g{ MC \j,N^ c \t^) = 9l (j - l,N^-3,t) + c , (2.10) 

where the simple-hairpin contribution to Co is found as 

c = 2p^ = 2/9 . (2.11) 

More complicated side branches will contribute also, but we observe (see Sect. 4) that 
reasonable changes of c in fitting to the data can be compensated by readjusting p . We 
thus by convention choose the value (2.11), which then fixes po to the value (2.5). We have 
checked that for microscopic times 10 < t( MC ) < 50, g[ MC ^ is well represented as 

g (MC) = 2 consf (t (MC) ) x > xpal/2 ^ 

9 

so that this choice of c is well justified. 

For the endsegments the correction is more important. These in a single move will jump 
a mean squared distance 

ci=2 , (2.12) 
and this motion is not taken into account in the theoretical model. We thus have 

g {MC)^ N (MC)^(MC) )=gi ^ N (MC)_^ t) + Ci (213) 

We have checked that the correction c\ — 2 precisely takes into account the difference in the 
motion of the end segment and the adjacent interior segment of the Monte-Carlo chain. 
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Though these corrections are microstructure effects, they cannot simply be ignored. In 
particular it is important to correct for the endsegment motion, since g[ M (1, N^ MC '\t^ MC '^) 
reaches values of the order 100 only for t^ MC > ~ 10 6 . 



Besides g[ MC ^ we also have measured the cubic invariant 

g[ MC \j,N^ MC \t^) 



E h*(t (MC) ) - r i>a (0) 



L \a=l 



1/2 



(2.14) 



It is easily checked that this function for an interior bead in our model reduces to the second 
moment of n(j, t) 



9i 



(j,N,t)=£ s (n%j,t) 



1/2 



(2.15) 



this relation holding as long as the bead stays in the initial tube. For the endsegment the 
expression is more complicated and given in appendix B of I . Again the relation among g\ 
and g[ M involves microstructure corrections, which, however, are more difficult to estimate 
and will not be considered. 



III. A FIRST INSPECTION OF THE SIMULATION RESULTS 
A. Comparison to previous work 

Reptation theory predicts power law behavior 

r t 1 / 4 , t < t < t 2 

9l (j,N,t) ~ {t/N) l '\ T 2 « t « T 3 (3.1) 
{ t/N 2 , T 3 <t . 

Here To = 0(N°) is the microscopic time, till the segment motion feels the constraining 
environment, T 2 = 0(N 2 ) is the equilibration time of the internal motion, and T 3 = 0(N 3 ) 
is the reptation time, needed for a complete destruction of the original tube. The last line 
of Eq. (3.1) identifies the diffusion constant of the chain as D ~ 1/N 2 . 

Evans and Edwards || introduced the above described Monte Carlo model to test these 
predictions. They used obstacle lattices of spacing m ■ £q, m < 10, and chains of length 
j\[( M c) < 30. The runs seem to extent up to t^ MC ^ ~ 10 3 . Clearly according to present day fa- 
cilities this is a fairly small scale simulation. Still within the scatter of the data the results for 
the smallest spacings m = 1, 2 seem to verify the predictions (3.1). In particular the authors 
observe D ~ iV 2 , T 3 ~ iV 3 , as well as t 1/4 -regimes and t 1/2 -regimes for the motion of the cen- 
tral segment. To check these results, Fig. 2 shows our data for g [ MC) ( N^ MC ^ /2, N^ MC \t^ MC ^>) 
in the common doubly-logarithmic representation. As is obvious, a t 1//4 -regime starts around 
f( M c) ~ xo 3 and barely is observable for N^ MC ^ = 80. It fully is developed only for larger 
chain lengths. A t 1 ' 2 -regime is not observable for N^ MC "> < 160. It may be present for larger 
chains but its unambiguos identification needs at least a further decade in time. Recall that 
these data are taken for the obstacle lattice of highest density, m = 1, as all our data. We 
conclude that the observation of amounts to a misinterpretation of the direct crossover 
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from the initial behavior, which roughly follows a t 1 / 3 -law, to free diffusion, as seen here for 
short chains. 

Deutsch and Madden || used the Evans-Edwards model with m = 1 to measure the 
diffusion coefficient D by following the center-of-mass motion of the chain: 

(R cm (t {MC) ) - R cm (0)) 2 ) t(M ^°° D tl ua > (3.2) 

Measuring chains up to length 100 they found D ~ (]y( MC )\ ; i.e., a considerably faster 

decrease than predicted by reptation theory. Our own data for the center-of-mass motion 

are shown in Fig. 3. We clearly can extract the diffusion coefficients up to iV (MC) = 80. For 
N {MC) > on i y 

an upper bound D < 0.2 (iV^^) 2 can be given. Our measured values 
of (jV (MC,) ) 2 • D are plotted against (jV (MC,) )~ 1/2 in Fig. 4. We also included data extracted 
from Fig. 3 of |J. Clearly the two sets of data are completely consistent. They nicely are 
fitted by the ansatz 

D = 0M(N( MC ^- 2 [l + 50(iV (MC) )- 1/2 l , (3.3) 



this form being motivated by Doi's work It, however, is clear that the range of chain 

lengths from 10 to 160 in Fig. 4 is insufficient to fully justify an ansatz leading to such 
a large first order correction. Still it shows that with chain lengths that presently can be 
reached, we are far from extracting the large- N limit of the diffusion coefficient. 



B. Time scales 

The reptation time T3 gives the time needed to destroy the tube. It may be defined in 
terms of the decay of the end-to-end vector correlations of the chain 0, but for the present 
work another definition is more convenient, both theoretically and experimentally. We define 
T 3 by the relation 

N 

L n max (T 3 ) = — , (3.4) 

which via Eq. (2.2) implies 

9l (0,N,T 3 ) = ((r (T 3 ) - r (0)) 2 ) 

= N = R 2 e , (3.5) 

where i? 2 is the mean squared-average end-to-end vector. Thus within time interval T 3 the 
endsegment has moved mean squared distance i? 2 . Our simulation data allow for the deter- 
mination of T 3 for N^ MC ^ < 160, and our experimental results together with the theoretical 
curve are shown in Fig. 5. Being based on an approximation the theory lies about 20 % 
above the datapoints. This suggests that our approximation underestimates n max {t). Both 
theory and data, however, consistently show that it needs chain lengths much larger than 
jy(MC) _ 200 to approach the asymptotic N 3 behavior. Indeed, the theoretical asymptote, 
calculated from Eq. (I 5.50), is found as 
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rf C) = 2.62(M MC )) 3 



(3.6) 



suggesting that chain lengths much larger than 10 3 are needed. This is quite consistent with 
our finding for the diffusion coefficient. 

Another important time scale of the model is the internal equilibration time T 2 of the 
chain. It gives the time a hairpin needs to diffuse over the whole chain, so that for t ^> T 2 
the motion of all beads is correlated. Theoretically T 2 can be identified with the Rouse time, 
i.e , the longest internal relaxation time of a free chain: pT 2 = N 2 /ir 2 . A practicable and 
precise experimental definition is not easy. Heuristically we could think of the time at which 
the motion of the central segment bends over from t xi ,xi xs 1/4 towards t X2 ,x 2 ^ 0.5, but 
this crossover is quite broad and the power law regimes are poorly defined for shorter chains. 
We thus here are content with the theoretical definition, valid for long chains: 

,T,= ( 7 ) (3.7) 

Transforming to the Monte Carlo time, we find 

T 2 (MC) = 1.66(iV( MC )) 2 , (3.8) 

where we ignore the difference among N, N^ MC ^ . Comparing to Figs. 2 or 3 we note that this 
marks the point, where definite deviations from the initial behavior can be seen. Combining 
Eqs. (3.6), (3.8) we find 

T (MC) 

!>=1.5*NW . (3.9) 

1 c 

Since both, J^ ) an d x^ MC \ mark fairly broad crossover regions, the ratio (3.9) must take 
values of order 10 3 before we can see the intermediate t 1//2 behavior. (Cf. the second line 
of Eq. 3.1.) Only the longest chain N = 640 shows a sufficiently large ratio ^ MC ^^ MC \ 
But then the total MC time 10 7 is not much larger than T^^ ^ = 0.7 • 10 6 . Even with the 
present data we thus have no chance to verify the t 1 / 2 -law. Qualitative, not quantitative, 
indications can be found however, as is discussed in sect. VI. These findings are completely 
consistent with typical results found in the literature for simulations of melts. 



IV. MOTION INSIDE THE TUBE 

A. Analysis of gi(^,N,t) 

In Eqs. (1.4.1), (1.4.2), (1.3.12) we have given the rigorous theoretical result for gi(j, N, t) 
for motion within the tube. To keep track of this condition we here denote this result as 
9i(j,N,t). 

9i (j, N, t) = (1 £ 2 p A(j, 0) V2 t 1 " *i( t))\ (4.1) 
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t 1 7V " 1 



,2 



irk 
2N 



- 2 (#C? + §)) (4 . 2) 



sm 2 ( ^ 



1 — exp [ —At sin 

\ ZiV / / J sin 2 , 2N , 

We use the variable i = pt, and we somewhat simplified the expression, the simplification 
being valid for pt > 1. The correction function Fi(j,t) reads (Eq. 1.3.26) 

It arises from the discreteness of the stochastic variable n(j,t) and is negligible for z > 25. 

As has been discussed in I, Sect. 4, for long chains and times so large that Fi(z) can be 
ignored, gi(j,N,t) takes the form (cf. Eq. I 4.10): 

g l (j,N,t) = (t*p y/HV%(^,^ ] j , (4.4) 

f 27T-^(^ po) l/2£l/4 |/7V2«1 



It is in this large time region that we determine the nonuniversal parameters. Specifically 
we get one relation from fitting the t 1//4 -plateau. 

^MOyl/Ag. ^ ^ ^ = 27r -3/4 (£ 2 po) l/2 T -l/4 j j<MO) ^ ^MC) ^ jjMC) ^ 

Since no t 1 / 2 -regime properly is reached by the data, we determine r, i^po separately by 
fitting to the crossover at t ~ T 2 , where the t 1//4 -regime terminates. In the fit we exclusively 
used data for the longest chain: N^ MC ^ = 640, so that we have a large region affected neither 
by initial effects nor by tube renewal. We find the values i^po = 1.23 and r = 6.092 • 10~ 2 
cited in Sect. II C. 

With £ 2 s p fixed, a variation of p only influences the argument of F\ in Eq. (4.1). In- 
creasing p we decrease the time range where iq is important. Since iq(0) = 1, we thus also 
increase the initial slope of gt(j,N,t). As mentioned in Sect. II C, these changes to some 
extent can be compensated by a change of Co that relates g[ MC ^ to g-y. Good fits can be 
obtained for 0.14 < c ^ 0.25, with varying p Q from 0.3 to 0.18. For t pa 10, gi(j, N, t) varies 
by about 15 %. Larger changes of the parameters lead to a mismatch in the curvature of 
the theoretical and experimental results. We finally fix all parameters by choosing c = 2/9, 
leading to p = 0.22. 

Fig. 6 shows our results for log 10 (y, N, tj / g ass (t)^ as function of log 10 (£), where 
9ass{t) gives the intermediate asymptotics defined by the first line of Eq. (4.5). This plot 
focusses on the region of t 1 / 4 -behavior. As mentioned above, the parameters have been fitted 
for M MC ) = 640, the remaining results involving no further fitting. This plot demonstrates 
the ability of our reptation model to explain the data. Deviations occuring for shorter chains 
and large time result from tube renewal, as is obvious from the values of T3 indicated by 
the ends of the full lines. Within a rough approximation such effects will be considered in 
Sect. V B. The discreteness correction Fi(z) is visible up to t ~ 10 3 , and for the shorter chains 
this initial range immediately crosses over to behavior dominated by tube renewal, with no 
indication of an intermediate t 1//4 -plateau. Furthermore, the t 1 / 2 -regime is not developed 
within the range of our simulations. 
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B. Analysis of <?i(f ,N,t) 

As has been stressed and illustrated earlier M, the cubic invariant <?i(y,iV, t) (Eq. 2.14) 
should show t 1/,4 -behavior even for very small time, with no discreteness correction. For 
motion inside the tube the theory yields 

MN,t) = (2e aPo A 1 (j,t)) 1 ' 2 , (4.7) 

with Ai(j, t) given in Eq. (4.2). Fig. 7 shows our numerical and analytical results, normalized 
to g aS s(t) = yf 9ass{t)- As expected, the t 1//4 -plateau is seen very clearly. Even for the 

shortest chain, N^ MC ^ = 20, there is an initial range of behavior close to t 1 ^. The plateau- 
value, however, systematically seems to lie 3-5 % below the theoretical value. We believe that 
this might indicate some small effect of the interaction among the hairpins, renormalizing 
the amplitudes but not changing the power laws. 



V. TUBE RENEWAL EFFECTS 

A. Motion of the endsegments 

Tube renewal cannot be treated rigorously, and our approximation, discussed in I, 
Sect. V, yields for the motion of the endsegment j — 

9l (0,N,t) = J2 -9*(0,N,s) (5.1) 
s =i s 

where ^(0, N, s) is given in Eq. (4.1). In Fig. 8 we have plotted <?i(0, N, t)/t l / A together with 
our data. Obviously our approximation reproduces the qualitative features of the data and 
performs not too bad on the quantitative level. It somewhat underestimates the mobility of 
the endsegment, consistent with our finding for T 3 (cf. Sect. Ill B), but the relative error 
decreases with increasing time. This is plausible since its origin lies in a mistreatment of 
the time-dependent correlations which decay on time scale T 2 . Both theory and experiment 
agree in exhibiting a very long initial transient, much longer than found for the central 
segment. Only the longest chains barely reach a t 1//4 -plateau. As mentioned in I, Sect. V C, 
this again is due to the discrete nature of the process. With £^p fixed, the theoretical 
curves again are fairly insensitive to po- We should recall that we have corrected the data 
by subtracting c\ — 2 from the motion of the endsegment (cf. Eq. (2.13)). Without that 
correction, the data in the initial range would be enhanced somewhat, starting around .54 
at t — 1. 

The theory predicts that gi(0, N, t)/gj(y, N, t) for T <^t <^T 2 reaches a plateau, which 
in our approximation is found at a value 4v^2- For t > T 2 this ratio should decrease again, 
since all segments start to move coherently. Fig. 9 shows our results for this ratio where for 
the central segment we in the theoretical result approximately took the tube renewal into 
account in the form discussed in the next subsection: , iV, i) = gi{ y , N, t) + 2g r (y , N, t). 
This figure should be quite characteristic for reptation. For t < T 2 it very clearly shows the 
enhanced mobility of the endsegment. Note that this ratio for the Rouse model would take 
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a value of 2 in the initial range and drop down to 1 for t ^> T 2 . Here it takes much larger 
values, coming close to our theoretical estimate. Clearly the asymptotic plateau value can 
not be estimated from the experiment. This would need much longer chains and longer 
times. For t > T 2 we see the expected decrease. For t > T 3 the ratio asymptotically should 
tend to 1. Also this is seen, but it is clear that it needs times t > 10 T 3 to approximately 
attain that limit. (In Fig. 9, T 3 for jV( MC ) < 160 can be taken from the endpoints of the 
theoretical curves.) 

We also measured the fourth moment #i(0, N, t). As we have discussed in appendix B of 
I, after some initial range the ratio «7i(0, N, £)/<7i(0, N, t) should take values 

gi(0,N,t) = f 1.0125 T <t<T 2 

gi(0,N,t) \ 1.085 T 2 <t<T 3 , 1 } 

where these results again are only approximate. Our data in all the regime T <C t C T 3 
scatter within bounds (1.025, 1.035), with no systematic trend observable. This is completely 
consistent with the estimate (5.2), if we take into account that within the limitations of our 
simulation the region T 2 <C t <C T 3 is not properly developed. 



B. Motion of an arbitrary segment 

The mean squared displacement of an arbitrary segment j not to close to the center of 
the chain for all t < T 3 can be written as 

9l (j, N, t) = 9i (j, N, t) + g r (j, N, t) , (5.3) 

where gi(j,N,t) is the contribution of motion inside the tube as given in Eq. (4.1), and 
g r (j,N,t) is the contribution of the tube renewal, which reaches segment j roughly at a 
time T R (j) defined by 

is n max (T R (j)) = j (5.4) 

(cf. I, Eq. 5.31). In I, Sect. V C and appendix C, we have evaluated g r (j,N,t) in a 
rough approximation, based on the distribution of n max (t) for a simple random walk, with 
hopping probability adjusted to our theoretical result for n max (t). Despite ignoring detailed 
correlation effects this approximation yields quite reasonable results as is illustrated in Fig. 
10. We there show the motion of segments j = 20, 40, 80 in a chain of length N^ MC ^ = 640. 
In particular for j = 80 we see the onset of the initial t 1 / 4 -plateau, which for i > 10 3 ' 5 is 
destroyed by the hairpins diffusing in from the nearest chain end (j = 0). This in itself 
would lead to another plateau at relative height \/2 (cf. I, Sect. IV), but before this can 
develop two further effects set in. Hairpins from the other chain end become important, 
too, (which implies t > T 2 ), and tube renewal effects are seen. Detailed inspection of Fig. 
10 shows small but significant deviations among experiment and theory, which must be due 
to our neglect of correlations in g r . 

Eq. (5.3) ceases to be valid close to the center of the chain, since it considers tube renewal 
only coming from one chain end. However, for t ^ T 3 processes where both chain ends within 
time interval t have made excursions of length w N/2 into the tube are rare. This suggests 
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to apply Eq. (5.3) also to the central segment. The result, which should underestimate the 
tube renewal effects, is given by the long dashed curves in Fig. 6. The short dashed curves 
follow by weighting g r (j, N, t) by a factor of 2, in trying to take into account the symmetry 
of the tube renewal process. The results look quite reasonable and clearly demonstrate that 
the deviations from <7i(y, N, t), as given by the full curves in Fig. 6, indeed are due to tube 
renewal. 

The motion of segments j = 80 for chains N^ MC "> = 160 and 640 is compared in Fig. 11. 
It for N( mc ^ = 160 shows the additional mobility due to a superposition of the effects from 
both ends. 



VI. QUANTITIES INVOLVING THE CENTER-OF-MASS 

Our theory at present has not been evaluated for quantities like the center-of-mass mo- 
tion, which involve two-bead correlations. We thus here only present our numerical results 
and compare to the qualitative predictions of reptation theory. 



A. Center-of-mass motion 

We have measured the correlation function 

g cm {t) = ((R cm (t) - R cm (0)) 2 ) , (6.1) 

where R cm (t) is the position of the center-of-mass of the chain. Reptation theory predicts 
the power laws 




Note that free diffusion sets in for t > T 2 , in contrast to the motion of the internal segment, 
where diffusional behavior is found only for t ^> T3 (cf. Eq. (3.1)). As shown in Fig. 3 our 
data reach the diffusional regime only for N^ MC ^ < 160, allowing for an extraction of the 
diffusion coefficients as discussed in sect. Ill A. Here we consider the short time regime. 

Fig. 12 shows the combination g cm (t) ■ N^^/i 1 / 2 . For N {MC ^ > 320 we indeed find a 
plateau. To the best of our knowledge this is the first time that the t 1//2 -behavior for the 
center-of-mass motion has been observed. The plateau seems to approach an asymptotic 
value close to 1.5, but the splitting of the curves even in the initial range indicates the 
existence of sizeable corrections to the iV-dependence. For shorter chains the large mobility 
of the chain ends ruins the t 1,/2 -behavior and a glance to Fig. 3 shows that effective power 
laws g cm ~ t x , I < x < .8, might be extracted, for t < T 2 . 

B. Motion of the central segment relative to the center-of-mass 

The correlation function g 2 (j,N,t) defined as 

g 2 (j,N,t) = (([ Tj (t) - R cm (t)] ~ M0) - R cm (0)] )^ (6.3) 
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measures the motion of segment j relative to the center-of-mass. For t <C T 3 the center-of- 
mass moves much slower than any specific segment, and thus 

g 2 (j,N,t)n 9i (j,N,t), t<T 3 . 

For i > T 3 , however, g 2 (j,N,t) saturates at some j-dependent value, that for j = N/2 
equals the mean squared radius of gyration R 2 , up to correction of order 1/N 2 . Reptation 
theory thus predicts 

N ( t 1/4 , T <t<T 2 

g 2 (-,N,t)~} (t/N) 1 / 2 , T 2 «t«T 3 (6.4) 
{ R] = N/6 , T 3 « t . 

Specifically in some intermediate range g 2 (^,N,t) increases with a larger effective power of 
t than in the initial range, and this phenomenon here is not mixed up with crossover towards 
free diffusion. Its observation thus is a clear signal of reptation. 

Fig. 13 shows our results for N^ IC ^ = 80, 160, 640, normalized to the t 1/4 -plateau g ass (t). 
The sequence of the three regimes (6.4) for iV( MC ) = 80, 160 is clearly seen, the data also 
saturating at R 2 , as expected. As also was to be expected, in the intermediate range 
T 2 < t < T 3 the power law g 2 ~ t 1 ! 2 is not fully attained, but still the data in this plot show 
a pronounced maximum. To our knowledge this is the first time that the intermediate t 1 ! 2 - 
regime has clearly been identified. In Fig. 13 we included theoretical curves for ^(y, N, t) 
(Eq. 4.1) to check whether g 2 (y, A?", t) for i < T 3 indeed equals (fc(~,iV, t). Taking into 
account that the data are not measured with very good statistics, being averaged over only 
10 independent runs each, we find a very satisfactory agreement. 



VII. CONCLUSIONS 

We have performed extensive simulations of the Evans-Edwards model up to chain 
lengths iV = 640 and 10 7 Monte Carlo time steps. Our simulation data show all features 
predicted by reptation theory, in particular: 

1) We find a strong increase in mobility of the endsegment, as compared to the central 
segment. 

2) The simulations of the motion g 2 of the central segment relative to the center-of-mass 
exhibit all three time regimes predicted by reptation, including the intermediate 't 1 / 2 '- 
regime. 

3) The crossover time T 2 to 't 1 / 2 '-behavior of g 2 coincides with the crossover time to free 
diffusion of the center of mass. 

These features clearly distinguish reptation from pure Rouse motion or a Rouse model with 
randomly spaced entropic traps [|H],|ll[ . 



4) We also found the celebrated t 1,/4 -law for motion of an inner segment, and the corre- 
sponding t 1//2 -law for the center of mass. However, we need chain lengths iV > 100 
and correspondingly long times to reliably identify such asymptotic laws. 
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5) Within the range of our simulations, the asymptotic iV-dependence predicted for the 
diffusion coefficient D, and the reptation time T 3 is not yet reached. From our results 
we estimate that chain lengths larger than N = 10 3 are needed to come close to 
asymptotics for these quantities. 

Our work shows that strong preasymptotic effects are an inherent feature of reptation. 
Such effects in fact dominate the chain length and time range covered by our simulations. 
The crossover regions here in particular cover all the region T2 < t < T3, masking in 
gi(j,N,t) the expected t 1//2 behavior for an inner segment. Also crossover from the initial 
behavior to the t 1//4 -law is very slow, a feature which we trace back to the discrete character 
of the basic dynamics. Indeed, this crossover is so slow that no t 1//4 -regime is seen for the 
endsegment. As Figs. 5-11 and 13 show, this is well explained by our theory. We indeed 
find very good agreement between our simulation data and our analytical evaluation of De 
Gennes' reptation model, also in the regions where the asymptotic predictions fail. To reach 
this agreement for the large variety of quantities and the large parameter range considered, 
we have adjusted the four parameters po£l,po,co and r = pt/t^ MC ' within the physically 
reasonable range. Our analytical predictions are exact, as long as tube renewal is negligi- 
ble. Deviations between theory and simulations are mainly due to our only approximate 
analytical evaluation of tube renewal. 

Will our numerical and analytical results be stable under a change of the microscopic 
structure of the system? This is a question of many different facets. It is expected that cer- 
tain types of time-independent disorder in the surrounding, entropic traps [T(| in particular, 



may ruin reptation all together. Also the consequences of relaxation of the surrounding, 
like in a melt, at present are not well known. Restricting ourselves to the original reptation 
scenario, i.e., to motion through an ordered array of obstacles, we should consider the effect 
of excluded volume interactions among the beads of the chain. In I, Sect. II C, we have given 
reasons why we believe this to be basically irrelevant here, changing only the embedding of 
the tube into real space as well as the time scale. More serious is the fact that both in our 
simulations and our analytical work we use a very narrow tube. Allowing for more degrees 
of freedom of the chain per unit spacing of the obstacle lattice, we certainly will increase 
microstructure corrections related to excursions of the bead considered from the center of 
the tube, and at the same time the discreteness corrections, playing such an important 
role in the initial time range within our model, will decrease. These effects to some extent 
may compensate each other. Indeed, in some preliminary simulation using a spring-and- 
bead chain in continuous space and obstacles of finite diameter in a regular lattice of wider 
spacing, we found results closely resembling the initial time range of our model presented 
here. (With this other model we, however, were unable to reach chain lengths and a time 
range where reptation predictions like the t 1//4 -law properly are found. Rather we stayed 
with effective t 1//3 -behavior familiar from previous work. With regard to the range of chain 
lengths and times this is completely consistent with the present findings.) We thus expect 
that certainly on the qualitative and presumably also on a semiquantitative level our results 
stay valid for wider tubes even in the initial time range. Of course a naive rescaling of our 
results in such a problem involving several time- and length-scales may be questionable. 

We finally should comment on consequences of this work for previous work on polymer 
motion through more realistic environments. Let us first consider polymer motion through 
a fixed disordered background of other chains, roughly modelling a gel. There the general 
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folklore tells us that the reptation scenario is valid. To examine this we appeal to ref. |12 
where the motion of a long chain (N = 200) was simulated. Fig. 6 of that work indeed 
shows more than one decade of t 1//4 -behavior for g 2 (y>-^£)> i.e , the motion of the central 
segment relative to the center-of-mass. However, the characteristic increase of the effective 
power at the end of the t 1//4 range, as shown in Fig. 13, is missing, this law fairly abruptly 
ending in saturation. Also g\ and g 2 deviate rapidly, g\ reaching only an effective t 1//3 -law. 
Both these observations are not compatible with our reptation results, and this conclusion 



is strengthened by a glance to the motion of the center-of-mass, as shown in Fig. 7, ref. |L2 
To our feeling this suggests that the published work shows disorder effects PH|JTT|JT3|] rather 
than reptation. 

With regard to polymer motion in melts the situation is even less clear. Again considering 
as example some extensive published work [FJJ, we note that many results shown there 



resemble our results found for short chains. In particular, the ratio <7i(0, N, £)/<7i(y, N, t) 
reaches values of order 2.7, i.e , larger than for a Rouse chain. For g 2 no tendency towards 
t 1,/2 -behavior is seen, but this might be due to the effective shortness of the chains, reaching 
only of the order of 6 entanglement lengths. Thus it is not unlikely that these results reflect 
'reptational' behavior of a very short effective chain. 

Clearly with regard to such more complicated systems much work still has to be done, and 
we hope to have contributed to this task by clearly exhibiting the quantitative consequences 
of the reptation model. 
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Figure captions 

Fig. 1 

An allowed hairpin- move and a forbidden kink-jump, illustrated for a chain embedded in a 
square lattice with obstacles in the centers of the cells. Our simulation uses the 3-dimensional 
version of this model. 

Fig. 2 

logio Ui MC) (t (MC) ,iv(^) )t (MC)^ as f unc tion of log 10 (t( A/c )) for N^ 1 ^ = 20,40,80,160,320,640. 
The straight lines correspond to power laws g\ ~ t 1 / 3 , t 1//4 , t 1//2 or g 1 ~ t. 

Fig. 3 

\og 10 R 2 cm ,R 2 cm = ^(R cm (t( CM )) -R cm (0)) 2 ^ as function of log 10 t^ MC \ Chain lengths 

N (MC) = 20,40,80, 160,320,640. The lines indicate the asymptotic behavior Dt {MC) (3.2). 
The thus determined D(N^ MC ^) is further evaluated in Fig. 4. 

Fig. 4 

( N {MC)y D ag f unct i on Q f (jyt^))" 1 / 2 , Dots: present work, ellipsoids: ref. g|. The straight 
line represents Eq. (3.3). 

Fig. 5 

(iV^) 3 T 3 (MC) as function of M MC ). The curve gives the theoretical prediction. The bro- 
ken line gives the theoretical asymptote. Points from our simulations. 

Fig. 6 

log 10 (#1 Cj,N,t \ I g a ss(t)^j as function of log 10 (f). Points are our data. From left: N^ MC ^ = 
20, 40, 80, 160, 320, 640. The full lines give the result (4.1), valid for motion inside the tube. 
The broken lines approximate tube renewal, as explained in sect. 5.2. All lines, except for 
iV > 320, end at T^(N). The horizontal broken line gives the intermediate asymptotics, here 
normalized to 1 by dividing through g ass = 27r~ 3 / 4 (i^Po) 1 ^ 2 t 1 ' 4 . 

Fig. 7 

As figure 6, but for the fourth moment g\ (y,-/V, t^j /g ass {t). Curves are for motion inside 
the tube. 

Fig. 8 

log 10 (#i(0, iV,*)/* 1 / 4 ) as function of log 10 (t). Data for N^ MC ^ = 20 - 640 (from left). The 
curves are calculated within the approximation (5.1) and end at T^(N). 

Fig. 9 

gi(0, N,t)/ gi (^Y,N,tj as function of log 10 (f). Data and theory like in Fig. 8. The broken 
line indicates the theoretical plateau value 4\/2. 

Fig. 10 

log 10 {gi{j,N,t)/g ass (t)) as function of log 10 (t) for j = 20,40,80, and N {MC ^ = 640. The 
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broken line gives the contribution gi(J,N,t), the full line is our result (5.3). 
Fig. 11 

\og 10 (gi(j,N,t)/g ass (t)) as function of log 10 (t) for j = 80 and M MC ) = 160,640. Curves 
give the full theory. 

Fig. 12 

gcm{t)N} MC) /i 1 ' 2 as f unc tion of log 10 (*) for = 20, 40, 80, 160, 320, 640 (from left). Val- 

ues of T 2 (N), calculated according to Eq. (3.7), are indicated by the arrows. 

Fig. 13 

log 10 [g 2 ^,N,tj j g a ss(tf) as function of log 10 (£). Chain lengths are indicated. The fat 
curves give (y, N, tj / g ass - The thin lines represent the asymptotic law g 2 = R 2 . The thin 
broken line illustrates a power-law g 2 ~ t 1 ! 2 . 
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